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We consider the random wetting transition on the Cayley tree, i.e. the problem of a directed 
polymer on the Cayley tree in the presence of random energies along the left-most bonds. In the 
pure case, there exists a first-order transition between a localized phase and a delocalized phase, 
with a correlation length exponent Vpure = 1- In the disordered case, we find that the transition 
0^ ' remains first-order, but that there exists two diverging length scales in the critical region : the 

, typical correlation length diverges with the exponent Utyp — 1, whereas the averaged correlation 

' length diverges with the bigger exponent Uav = 2 and governs the finite-size scaling properties. 

04 ' We describe the relations with previously studied models that are governed by the same "Infinite 

(-H , Disorder Fixed Point". For the present model, where the order parameter is the contact density 

■ 6l = la/ L (defined as the ratio of the number la of contacts over the total length L ), the notion 
of "infinite disorder fixed point" means that the thermal fiuctuations oi 9l within a given sample, 

, become negligeable at large scale with respect to sample-to-sample fluctuations. We characterize 

the statistics over the samples of the free-energy and of the contact density. In particular, exactly 
at criticality, we obtain that the contact density is not self-averaging but remains distributed over 
the samples in the thermodynamic limit, with the distribution VtcX^) ~ l/(7r-^S(l — 9)). 

C/3 . I. INTRODUCTION 

In pure phase transitions, the approach to criticality is usually governed by a single correlation length exponent 

^ ' V that describes all finite-size scaling properties. In the presence of frozen disorder however, the distribution of 

^ ^ correlation functions may become very broad and disorder-averaged values can become atypical and dominated by rare 

' . ' events. This phenomenon has been first studied in one dimensional classical spin systems [l|, SjlH, where correlation 

^ functions can be expressed as product of random numbers, and have been then found in higher dimensional systems, 

Q such as the 2D random g-state Potts model [5] and the random transverse field Ising chain [4] (which is equivalent 

Q \ to the classical 2D McCoy- Wu model) . It is important to stress that both the typical and averaged correlations are 
actually important, depending on the physical quantities one wants to study [l|,H, It turns out that close to a 

I ' phase transition, the typical correlation length ^typ and the disorder-averaged correlation length £^av may have different 

^ . critical behaviors. The best understood example is the random transverse field Ising chain, which has been studied 

f> — ' in great detail by D.S. Fisher via a strong disorder renormalization approach jj] to obtain utyp — 1 and Vav = 2. 

0^ These two exponents vtyp — 1 and Vav = 2 also occur in other disordered models that are described by the same 

■ "Infinite Disorder Fixed Point" (see the review [S] and references therein). In the present paper, we present still 
' another realization of this "Infinite Disorder Fixed Point" with vtyp = 1 and Vav = 2, as a random wetting transition 

y—^ on the Cayley tree. 

\_ Our physical motivation to consider such a model of random wetting on the Cayley tree was to better understand 
0^ . the similarities and differences with two other types of models involving directed polymers and frozen disorder : 

' (i) for the problem of the directed polymer in a random medium on the Cayley tree 0, there exists a freezing 
^ , transition towards a low-temperature phase of finite entropy, where the polymer is essentially frozen along the random 
optimal path. Finite-size properties in the critical region [3, have revealed the presence of two distinct correlation 
^ . length exponents v = 2 and v' = \. In the random wetting model considered in the present paper, the difference is 

■ that the the random energies are not on all bonds of the Cayley tree, but only on the bonds of the left-most path, so 
that the phase transition corresponds to a freezing along this boundary path. 

(ii) for the wetting [To'| and Poland-Scheraga model of DNA denaturation with a loop exponent c > 2, where 
the corresponding pure transition is first-order, we have found numerically [l2l ITsj that in the presence of frozen 
disorder, the transition remains first order, but that two correlation length exponents vtyp = 1 and v — 2 appear. The 
random wetting model considered in the present paper corres pon ds to the limit of loop exponent c — > 00 (since loops 
do not exist on the Cayley tree) of the model studied in Ref. [12| with the boundary conditions bound-unbound (see 
more details on the model in section 2.1 of Ref. p^). 

The paper is organized as follows. In Section|TIl we introduce the model and the interesting observables. In Section 
mil we discuss the statistical properties of the free-energy and conclude that the wetting transition remains first-order 
with a jump of the energy density. In section IIVI we study the statistics of the attached and detached lengths, and 
give the consequences for the statistics of the order parameter in the critical region. We summarize our conclusions in 
IVII Appendix A contains for comparison the analysis of the properties of the wetting transition in the pure case. In 
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Appendix B, we discuss the transition temperatures of the moments of the partition function in the disordered case. 

II. MODEL AND OBSERVABLES 
A. Definition of tlie model 




FIG. 1: In this paper, we consider a Cayley tree of branching ratio K and a number L of generations (K = 2 and L = 4 on 
the Figure), and we study the statistical physics of a directed polymer starting at the root in the presence of random energies 
(—Ml, —112, .., —ul) along the left-most bonds. A configuration of the directed polymer is characterized by an attached length 
Q < la < L along these left-most bonds and a detached length Id = L — la {la = I and Id ~ S on the Figure). The delocalized 
phase corresponds to a finite la as L ^ +oo, whereas the localized phase corresponds to a finite Id as L ^ +oo. The order 
parameter is the contact density 6 = la/L. 

We consider a directed polymer of length L on a Cayley tree of branching ratio K. The total number of directed 
walks is simply 

Yl - (1) 

The corresponding delocalized free-energy reads 

pdeioc = ^tIuYl = -TL In K (2) 

and corresponds to an entropy of (In K) per bond. 

We now consider the random wetting model, where the left-most bonds are characterized by independent random 
energies {—Ui) with i — 1, 2, .., L (see Figure[T]). As an example, they can be drawn from the Gaussian distribution of 
average value uq > and variance 

p{u) = , e 



(3) 



B. Recurrence satisfied by tlie partition function 



To write a recurrence for the partition function, it is convenient to label the random energies from the bottom, 
instead of the labelling from the root shown on Fig. [1] 

Ui = UL-i (4) 

The partition function then satisfies the simple recurrence 

ZL+i=e^''^^'ZL + {K-l)YL (5) 
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with the initial condition Zq = 1. It is convenient to introduce the ratio with respect to the free partition function of 
Eq. □ 



The recurrence then becomes 



with the initial condition i?o = 1. 



Rl^i = '-^Rl + (7) 



C. Free-energy difference with the delocalized state 

The difference between the free-energy Fl = —TlnZ^ and the delocalized free-energy of Eq. [2] is given by 

pd^ff ^p^_ pdeioc ^ _j, Zl+TIiiYl = -T In Rl (8) 
The excess free-energy per monomer due to the wall in the thermodynamic limit L — > +oo 

_ pdeioc 1^^^ 

foo = hm = -T hm —— (9) 

L — >+oc Li L — >+oc Li 

characterizes the wetting transition : it vanishes in the delocalized phase foo{T > Tc) ~ and remains finite in the 
localized phase foo{T < T^) < 0. 

D. Probability distributions of the attached length la and of the detached length 

Since there is no loop on the Cayley tree, the directed polymer cannot return to the wall after leaving it. A 
configuration can be thus decomposed into a length la attached to the wall starting at the root and a length l^ = L — la 
detached from the wall (see Fig. [1]). The thermal probability QL{la) to have exactly la links attached to the wall in 
a given disordered sample reads 



Ui 



The numerator is the contribution to the partition function of the configurations that leave the wall after exactly la 
steps. The denominator is the partition function 

la 

Zl=Y.^ »=i [5ia,L + {K- l)K'^-'^-^0{la < L)] (11) 
so that the normalization of the probability distribution QL{la) reads 

L 



J2 = 1 (12) 



la=0 

It will be convenient to introduce also the probability distribution of the detached length Id ^ L — la 



Ui 



4 



E. Order parameter : contact density 

A convenient order parameter of the transition is the contact density, which is directly related to the attached 
length la introduced in Eq. (TU] 



(14) 



In particular, its thermal average in a given disordered sample is determined by the first moment of the distribution 
QLiL) of Eq.m 



< >= = y laQdla) 



(15) 



la=0 



In the thermodynamic limit L — > oo, it is expected to remain finite in the localized phase and to vanish in the 
dclocalizcd phase. 

F. Internal Energy 

The internal energy 

EL = -dp In Zl = -dp In Rl (16) 
satisfies the recursion (see Eq. [7] for the corresponding recursion for the reduced partition function Rl) 

^f'^^ + ^RLi-UL + l + EL) 



El+1 



(17) 



and the initial condition Eq = 0. In terms of the attached length la introduced above (see Eq. [TO)) . it can also be 
written as 



L /la \ P^Ui 

E E "0 -1 [5ia,L + {K- l)K'^-'--'dila < L) 



El = 



Zl 



The energy per monomer in the thermodynamic limit L +oo 



Boo = hm — 

L^ + oo L 



(18) 



(19) 



remains finite in the localized phase e^oiT < Tc) > and vanishes in the delocalized phase eoo{T > Tc) — 0. Its value 
can be obtained via the usual thermodynamic relation from the excess free-energy f^o introduced in Eq. [5] 



(20) 



III. STATISTICAL PROPERTIES OF THE FREE-ENERGY 

In this section, we discuss the statistics of the free-energy difference of Eq. [8] which depends only on the variable 
Rl introduced in[6l We first explain that Rl is a Kesten random variable and describe its statistical properties. 



A. Analysis of the variable Rl as a Kesten random variable 

The recursion of Eq. [7| takes the form of the recurrence of Kesten variable [l3| 

Rl+1 — a,L+iRL + b 



(21) 
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where b = {K — 1)/K is a. constant and where the 

a. = (22) 

are independent identically distributed random variables. The specific structure of a Kesten variable Rl consists in 
a sum of products of random variables 

L / L L 

RL^Y[a, + bil + Y,l[^^] (23) 

1=1 \ j=2 t=j 

This discrete form has for continuous analog the exponential functional 

TZl = / dxe^^^ '^y^^y^ (24) 
Jo 

where {F{xy\ is the random process corresponding to the random variables (Ina;) in the continuous limit. 

This type of random variables, either in the discrete or continuous forrns, appears in a variety of disordered systems, 
in particular in random walks in random media [isl . [l6l [itI . in the classical random field 

Ising chain [2^, [23| , and in the quantum random transverse field Ising chain 28 , 291 30 | . Many properties have been 
thus already studied in great detail in these previous works. In the following, we cite and translate these known 
results for the present context. 

To describe the universal results near the critical point for large samples, there are only two relevant parameters : 

(i) the first important parameter is the averaged value 



Fa = \na^^ l3u-\nK (25) 



The critical point corresponds to the vanishing condition Inoj = 0, so that the critical temperature of the present 
random wetting model reads 

(In particular, for the Gaussian distribution of Eq. [3l the critical temperature is Tc — UQ/{\nK) and is independent 
of the variance A^.) 

(ii) the second important parameter is the variance of (Ina^) that will be denoted by 

2a = (InaO^ - (h^)' = /? - {ulf) (27) 

N 

to keep the same notations as in Refs. [l^, [2^, [s^. In particular, from the Central Theorem, the sum S^v = y~^(ln aj) 

i=l 

is then distributed asymptotically for large N with the Gaussian distribution 

Pn{J:n) ^ e- (28) 



B. Delocalized phase T > 

For T > Tc, the parameter of Eq. [^S] is negative Fq < and the random variable Rl remains a finite random 
variable as L ^ 00. Its probability distribution Poc {Roo) is known to display the power-law decay [13] 

/'oo(i?oo) „ (X (29) 

where the exponent /i(r) > is determined by the condition 

i = «r = i^ (30) 
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In conclusion, the difference of Eq. [8] between the free-energy and the delocahzed free-energy F'lf^°'^ remains a 
finite random variable as i ^ cxd 



pd^f^ = p _ p. 



delo 



pMff ^ R 



and its distribution presents the following exponential tail for Fdiff —f — oo (sec Eq. 



Fd.ff- 



Example : for the case of the Gaussian distribution of Eq. [31 the condition of Eq. [30] becomes 



1 = 



and the exponent fi reads 



2T{T\nK -uZ) _ 2{ln K)T{T - T^) 
A2 ^ A2 



(31) 



(32) 



(33) 



(34) 



C. Critical point T = 

Exactly at criticality T = Tc, the parameter of Eq. [25l vanishes Fq = 0. It is then known (sec for instance [2^.[23l[3C 
and references therein) that the leading behavior is of order 



lni?L(rc) ~ 2{aLY/^Wc 

L— *oo 



(35) 



where Wc is a positive random variable of order 0(1). Moreover, it can be shown (see [30l | and references therein) that 
the distribution of Wc is the half-Gaussian distribution (cf Eqs (24) and (32) in [30|) 



P{wc)^e{w, > 0)^e- 



(36) 



D. Localized phase T < Tc 

For T < Tc, the parameter of Eq. [25] is positive _Fo > and the random variable will grow exponentially in L 

~ (37) 



\nRL{T <Tc) ~ LFo + 2VaLw 

L — >oo 

The random variable w of order 0(1) is distributed with the Gaussian distribution (cf Eqs (24) and (31) in [30j ) 

1 ,„2 



P{w) 



(38) 



E. Summary of the critical behavior of the free-energy and energy 

difference 

the following critical behavior 



In conclusion, the difference of Eq. [5] between the free-energy Fl and the delocalized free-energy Ff^'^^°'^ displays 



F\ 



f"{T > Tc) 



= Fdiff of order 0(1) 



Ft"^^{Tc) = -Tc2{aL) 



1/2, 



with > of order 0(1) 



pdiff^j, ^ ^ _^- _ y _ T2VaLw with w of order 0(1) 



(39) 
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where the distributions of the random variables Fdi ff,Wc and w have been given above in Eqs [32l [36] and [38l In the 
thermodynamic Umit L +00, the excess free-energy per monomer due to the wall defined in Eq. [5] is thus given by 

/oo(r>Tc) =0 (40) 

foo{T<Tc) =-{u-T\nK)^-{T,-T)\nK (41) 

The thermodynamic relation of Eq. [201 then gives the behavior of the energy per monomer 

eoc(r > T,) =0 (42) 
eoo (T < T,) ^ -T, In K = -u (43) 

The transition is thus first-order, with a jump of the energy per monomer. 

IV. STATISTICAL PROPERTIES OF THE ATTACHED AND DETACHED LENGTHS 

In this section, we discuss the statistical properties of the attached length la or of the complementary detached 
length Id whose distributions have been introduced in Eq. [TUj and Eq. [T^j respectively. This allows to analyse also the 
statistics of the contact density of Eq. [141 which represents the order parameter and of the internal energy of Eq. [T8l 

A. Typical behavior of QL{la) in the delocahzed phase 

In the delocalized phase, the attached length la will remain a finite random variable as L — > 00. The typical decay 
of Qhila) is then given by (Eq. [TU)) 

\nQL{la) ^ piaU+iL -la)\nK + 1x1-^—^-1^ (44) 

K 

~-la{\nK-pu)+\n^-^^—^ + {L\nK-WZL) (45) 

K 

i.e. the decay is governed by the typical correlation length 

-^=\nK-(M (46) 
that diverges with the correlation length exponent 

Vtyp = 1 (47) 

as in the pure case ( Vpur = 1 of Eq. IA8|) . 

B. Typical decay of PL{ld) in the localized phase 

In the localized phase, the detached length l^ will remain a finite random variable as L ^ +00. The typical decay 
of PL{ld) is then given by 

In PL{ld) ^ (3(1 - ld)u + ld\nK + In ^^—^-hTZ^ (48) 

K 

c:^ -ld{f3u- In K) + \n^^—^ + {L(M-h[z2) (49) 

K 

i.e. the decay is governed by the typical correlation length 

-^=/3u-lnfC (50) 
that diverges with the correlation length exponent 

Vtyp = 1 (51) 

as in the pure case ( i/pur = 1 of Eq. IA17P . 
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C. Statistics over the samples of the probability = Qhila ~ 0) of zero contacts 



To understand why typical and disordered averaged behaviors can be different for the distributions Qhila) and 
Phi^d)^ it is convenient to consider first the particular case of the probability of zero contacts 

U.QaL-»)-Pai.-L)^^Ji^^JJi^±- (52, 

because it depends only on the Kesten variable Rl whose statistics has been discussed in detail in section Hill We 
also note that ql coincides with the expression of the flux in the Sinai model, and we refer to where many 
exact results have been derived for its probability distribution (see in particular Eq. (6.10) in [l^l)- In the following, 
we concentrate on the critical region. From the statistics of Rl described in section [Till one has in particular that 

(i) in the delocalized phase, qL remains a finite random variable as L — > +oo 

(ii) in the localized phase, the probability qL becomes exponentially small in L 

(ii) exactly at criticality, the variable qL behaves as qL ~ g-wc2\/aL -^j^gj-g jg g, positive random variable of order 
0(1) distributed with the half-Gaussian distribution of Eq. [Mj As discussed in [s^l, an important consequence of the 
critical statistics of Eql36| is that the averaged value is dominated by the rare samples having an anomalously small 
value of the scaling variable < w < and since the probability distribution P(wc) is finite at Wc = 0, one 

obtains the slow power-law decay 



gL(re) (X -j= (53) 

whereas the typical behavior is given by the exponential decay 

In gL(Tc) cx -VZ (54) 

L — >oo 

We refer to [30| and references therein for the discussion of similar examples of averaged values dominated by rare 
events in the field of quantum spin chains. 

D. Notion of " Infinite Disorder Fixed Point" 

The probability distribution QL{la) of the attached length la in a given sample has the form of a Boltzmann measure 

Qdla) = (55) 



E 



over the random walk potential 

la. la 

i=l i=l 

The continuum analog of Eq. [55] is then the Boltzmann measure over a Brownian potential (see Eq. 10 in ^2^). We 
refer to (25j and references therein for more detailed discussions, but the essential property is that the Boltzmann 
measure over a Brownian valley is very concentrated (in a region of order 0(1)) around the absolute minimum of 
the Brownian potential. Within our present notations, this means that in a given disordered sample at a given 
temperature, the thermal probability distribution Ql^o) is concentrated in a finite region around the point where 
the random potential of Eq. [51] reaches its minimum on the interval < < i 

QL{la)^5{la-l*a) (57) 

This means that in a given sample, the contact density of Eq. [TJ] is essentially given by its thermal average 



Ol ^< Bl >= I 



(58) 
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E. Sample-to-sample fluctuations of the contact density 

We now discuss the sample-to-sample fluctuations of Z* and of the contact density. The disorder average of the 
probabihty distribution QL{la) can be obtained as 



QLila)^6ila~l*)^nL{la) (59) 

where TTLi^a) the probabihty distribution of the minimum I* over the samples. 

Taking into account the two parameters Fo (Eq. [25|) and a (Eq. [27]) that characterize the large scale properties of 
the random walk U{la), the probability distribution 7ri(/*) of the minimum /* can be obtained as 

■^L(L) = —f (60) 

/o^rfZV(^o)(L-?*)V.(-^o)(/*) 



where 



= / dUG^^"\U,l) (61) 



and where G'''^°\U,1) represents the probability for a biased random walk to go from (0,0) to {U,l) in the presence 
of an absorbing wall at [/ = 0~ (see for instance j3(| for very similar calculations) 

The critical behavior can be analyzed as follows. 

1. Sample-to-sample fluctuations at criticality 

At criticality where Fq — 0, the function of Eq. [61] reads 

r+°° jj 2 1 

rP^'>=%l) ^ / dU ,,3/2 ^'^ = ^ (63) 

The distribution of Eq. [60] for the position I* then reads 

QUI*)-^^l'\i:)= /J^ (64) 
n^l*{L-ll) 

As a consequence, we find that the contact density 6* ~ is distributed at criticality over the samples with the 
law 

VtAS) = , ^ (65) 

2. Sample-to-sam.ple fluctuations in the delocalized phase 

In the delocalized phase where Fq < 0, the position I* of the minimum on [0, L] of the Brownian potential of positive 
drift remains finite as L — > +cx), and the probability distribution of Eq. [BO] can be obtained as 



'{Q^L^oo (66) 

with 



,<-.,(U . r,u-^^y^ - r (67, 
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In particular, it decays exponentially as 



with the correlation length 



Qi^a) =^f>^^^(/a) oc e-^'" (68) 
4cr 1 

^""Tf T^T+ {T-T,y ^^^^ 

This should be compared with the typical decay of Eq. US] As a consequence, the finite-size scaling properties in the 
critical region involves the correlation length exponent 

lyps = 2 (70) 
and not the typical correlation exponent that appears in Eq. HTl 

3. Sample-to-sample fluctuations in the localized phase 

In the localized phase, where Fq > 0, it is the detached length l'^ = L — I* that remains finite as L ^ +00, and the 
disorder-averaged probability of the detached length Id can be obtained as 



PM = <^^^(/: = L-Q (71) 

dllp(^o)(ld) 



with 







In particular, it decays exponentially as 



PL{ld) oc e-^^^ (73) 

id— »+oo 

This should be compared with the typical decay of Eq. 1491 Again, this means that the finite-size scaling properties 
in the critical region are governed by the correlation length exponent vps = 2 of Eq. [7D1 

V. DISCUSSION 
A. Lack of self-averaging of the order parameter at criticality 

Outside criticality, the contact density is self-averaging in the thermodynamic limit L -\-oo 

VtktAO) =5(0-1) (75) 

Vt>tAO) (76) 

but exactly at criticality, we have obtained that the contact density 6 remains distributed at criticality over the 
samples with the law of Eq. [65] 

VtAS) = , ^ (77) 

The lack of self-averaging of densities of extensive thermodynamic observables at random critical points whenever 
disorder is relevant has been studied in 31, 32, 33]. The main idea is that off-criticality, self-averaging is ensured by 
the finiteness of the correlation function ^, that allows to divide a big sample L ^ £,{T) into a large number of nearly 
independent sub-samples. However at criticality where the correlation length diverges ^ — -\-oo, even big samples 
cannot be divided into nearly independent sub-samples, and one obtains a lack of self-averaging. The cases considered 
in [m [12, [13 were second-order phase transitions, but the result of Eq. [51] means that the contact density is not 
self-averaging at the random first-order transition studied in this paper. 



with the correlation length 
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B. Physical interpretation of two correlation length exponents 

The presence of two different correlation length exponents Vav — 2 and i^typ = 1 whenever the disorder is equivalent 
to a Brownian potential U{x) is well known in other contexts, in particular for random walks in random media [T9l.[25j. 
for the quantum transverse field Ising chain 0, [s^l and more generally in other models described by the same "Infinite 
Disorder Fixed Point" (see the review Q and references therein). The physical interpretation is the following [3| : 

(i) the first length scale corresponds to the length £^typ where the mean value [U{L) — U{0)] — FqL is of order one, 
which yields 

6yp ~ with l^typ = 1 (78) 

(ii) the second length scale corresponds to the length where most of the samples indeed have {U{L) — f/(0)) ~ 
FqL ± V aL of the same sign of the mean value, i.e. the scale V aL of the fluctuations should be of the same order of 
the mean value, which yields 

£,av ^ with Vav = 2 (79) 

We have found that the finite-size scaling properties for disorder-avera ged values over the samples are governed by 
i^FS = ^au = 2 which actually saturates the general bound vps ^ '^/ddis l2J]i where ddis is the dimensionality of the 
disorder (here ddis = l)i whereas the typical exponent is smaller Vfyp = 1. This possibility of a first order transition 
that remains first order in the presence of quenched disorder has been already discussed in [sj] and in Sec. VII A of 
Ref. ^. 



VI. CONCLUSIONS AND PERSPECTIVES 



In this paper, we have studied the random wetting transition on the Cayley tree. We have obtained that the 
transition which is first order in the pure case, remains first-order in the presence of disorder, but that there exists 
two diverging length scales in the critical region : the typical correlation length diverges with the exponent vtyp — 1, 
whereas the averaged correlation length diverges with the bigger exponent Vav = 2 which governs the finite-size scaling 
properties over the samples. We have describe the relations with the previously studied models that are governed by 
the same "infinite disorder fixed point" . We have given detailed results on the statistics of the free-energy, and on the 
statistics of the contact density 9 = la/ L which constitutes the order parameter in wetting transitions. In particular, 
we have obtained that at criticality, the contact density is not self-a veraging but remains distributed over the samples 
in the thermodynamic limit, with the distribution VtS^) — l/(^\/ ^(1 ~ ^))- 

As explained at the beginning, our physical motivation to consider such a model of random wetting on the Cayley 
tree was to better understand the similarities and differences with two other types of models involving directed 
polymers and frozen disorder described in the points (i) and (ii) of the Introduction. The results given in the present 
paper suggest the following conclusions : 

(i) The freezing transition of the directed polymer in a random medium on the Cayley tree [7] is of a very different 
nature, since it is the exponent v = 2 which governs the free-energy singular part (instead of vtyv ~ 1 here). The 
appearance of a smaller exponent v' = 1 in some finite-size properties in the critical region [^, |0| should be then 
explained with another mechanism. 

(ii) Let us now compare with the wetting and Poland-Scheraga model of DNA denaturation with a loop exponent 
c > 2 that we have studied numerically in [l^, Ell • We find that the random wetting model considered in the present 
paper ( that corresponds to the limit of loop exponent c — > oo of the model studied in Ref. fT5| ) presents very similar 
critical behavior : in both cases, the transition remains first order with Vtyp = 1 and a finite contact density, but the 
finite-size scaling properties over the samples involve the exponent vps — 2. We thus hope that the exact results 
obtained here for the special case c — > oo will help to better characterize the transition in the region of finite-loop 
exponent 2 < c < -l-oo. 



APPENDIX A: PROPERTIES OF PURE WETTING TRANSITION m = uo 



In this Appendix, we describe the critical properties of the pure wetting transition on the Cayley tree, to compare 
with the disordered case considered in the text. 
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1. Finite-size partition function 

When all the energies Ui take the same value uq, the ratio of Eq. [6] simply reads 

L 



K K 



1 



/3"o\^ {K-\) {^) -1 



K ) \ K ) \ K 

L 



K I K (e^ _ I 



rypure 



X ^. (Al) 

K 

The critical temperature Tc then corresponds to the point where the energy per link uq exactly compensates the 
entropy per link (InK) of the free walk 

g/3eUo -> TP"™ = (A2) 



2. Delocalized phase e''"" < K 

In the delocalized phase, the leading term of Eq. I All remains finite 

Rl-iT>T.)^^^^^^^-,... (A3) 

and the excess free-energy per monomer due to the wall in the thermodynamic limit L +oo (Eq. [5]) vanishes 

/oo(T > T,) = (A4) 
The finite-size behavior of the full free-energy difference of Eq. [8] reads near criticality 

pdif f(T>T^^FLiT>T^~ Ff"'"" = -TIiiRl = -T In , ''^ ^ P , ~ -Tc In — ^ — ( A5) 

In the delocalized phase, the probability distribution QL{la) of the attached length la of Eq. [10] remains finite as 

L ^ oo 

QL{la)=PL{h = L-la) - l___ (Xe-/5"o) ^=Q^{Q (A6) 

L~too \ A / 



It decays exponentially 



where the correlation length reads 



QooiL) ^ e i^^io.iT) (A7) 

la — >- + 00 



e:iMT > = ^^^^ iT~Tcr.^r. -^^^ -p«- - 1 (^s) 

The thermally averaged contact density < 6*^ > of Eq. [15] vanishes as 

1 °° 1 
< > ,-00 Z E l^Q-(^^) - L(i^e-^-o - 1) (A9) 

i.e. in the critical region, it behaves as 



<'''>.2.^I(5:3« (AID) 
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3. Localized phase e''"" > K 



In the localized phase T < T^"""^, the leading term of the partition function of Eq. lAll is exponentially large in L 



1 



0(1) 



L^oo \ K 

The excess free-energy per monomer due to the wall in the thermodynamic limit L —>■ +oo (Eq. [5]) reads 

foo{T < T,) - -{uo - Tlnif) = ^^TPure __T)\nK 



(All) 



(A12) 



i.e. it vanishes linearly as the temperature T approaches the critical value J'P'"''^. The transition is thus first order. 
The internal energy of Eq. [11] becomes from Eq. lAllI 



^ ' L^+oc {ef^o - K) (e/5«o - 1) 

The intensive energy of Eq. [TH] remains constant in the whole low-temperature phase 



(A13) 



(A14) 



and presents a jump at criticality. 

The detached length Id remains a finite random variable in the thermodynamic limit L — > cx3 with the following 
distribution ( Eq. [T3)) 



Pooild 

It decays exponentially 

where the correlation length diverges as 



(K-l) 



P^{ld) ~ e ^ 



— tR^ n \ ^ 7^ \^ ^^^^ '^P"'-'^ ^ ^ 



(A15) 



(A16) 



(A17) 



4. Critical point e*^-"" = K 

Exactly at criticality, the ratio of Eq. lAll reads 



so that the free-energy difference is logarithmic in L 

The recurrence of Eq. [T7] for the energy becomes 

KRl 



(K-l) ■ 
1 + ^L 



L — >oo 



-Tr In L 



El+1 — 



KRl + iK-l) 



K+{K-l)L , 

\—Un + El A = : TT ziun + ErJ 



leading to 



K +{K - 1)L 



K 

{K-l)L 



(A18) 



(A19) 



(A20) 



(A21) 
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so that the energy per monomer in the thermodynamic Hmit L — > cx3 is 



2 



(A22) 



i.e. exactly in the middle of the jump between e^'^^{T < Tc) 
The probability distribution Qhila) of the attached length 



-Mo (Eq. EHI and eP^'-'^iT > = 0. 

for a system of size L becomes at criticality 



Qdla) 



+ ^^eila < L) 



1 + 



K 



L 



(A23) 



i.e. it is essentially flat over all values < la < L. The two phases coexist with proportions [lajld = L — la). 
So the value of Eq. IA21I simply means that the domain wall between the two phases is on average at the middle 
^ Id ^ but the domain wall is actually anywhere on the interval < < L with a flat distribution. 



APPENDIX B: DISORDERED CASEj_TRANSITION TEMPERATURES r„ FOR THE AVERAGED 

MOMENTS Zl OF THE PARTITION FUNCTION 

Since in other disordered models, one sometimes consider the series T„ of critical temperatures for the moments Z'^ 
of the partition function, we discuss their properties in the present Appendix. For n — \, one sees that the annealed 
partition function satisfies the recursion (see Eq. [7|) 

{K - 1) 



i?L+i = -j^Rl + (Bl) 

This is equivalent to the pure case recursion with the change e'^"" e''" so that the annealed critical temperature is 
determined by the condition (see Eq. IA2[) 



6*31" = K (B2) 

i.e. in terms of the exponent /i(r) introduced in Eq. I30| this corresponds to the condition 

m(Ti) = 1 (B3) 

More generally from the power-law of Eq. I29[ it is clear that the transition temperature r„ for the moments Z" of 
order n is determined by the condition 

M(Tn) = n (B4) 
For instance, for the case of the Gaussian distribution of Eq. [H the expression of Eq. [M] yields 



1 + ,/l + 271^41^ 



Tn - Tc ^ (B5) 
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